Dynamical mechanism of anticipating synchronization in excitable systems 
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We analyze the phenomenon of anticipating synchronization of two excitable systems with unidirectional 
delayed coupling which are subject to the same external forcing. We demonstrate for different paradigms of 
excitable system that, due to the coupling, the excitability threshold for the slave system is always lower than 
that for the master. As a consequence the two systems respond to a common external forcing with different 
response times. This allows to explain in a simple way the mechanism behind the phenomenon of anticipating 
synchronization. 



The synchronization of nonlinear dynamical systems is a 
phenomenon common to many fields of science ranging from 
biology to physics IB, and it has been an active research sub- 
ject since the work by Huygens in 1665. Recently, the syn- 
chronization of chaotic systems in a unidirectional coupling 
configuration has attracted a great interest due to its potential 
applications to secure communication systems yy. Particu- 
lar attention has been payed to the so-called anticipating syn- 
chronization regime, an idea first proposed by Voss in |3J. He 
showed that, in some parameter regions, two identical chaotic 
systems can be synchronized by unidirectional delayed cou- 
pling in such a manner that the "slave" (the system with cou- 
pling) anticipates the "master" (the one without coupling). 
More specifically, the coupling scheme proposed in 3 for the 
dynamics of the master, x(t), and slave, y{t) is the following: 
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where y T = y(t — r). For appropriate values of the delay 
time t and coupling strength K, the basic result is that y{t) « 
x(t+r), i.e. the slave "anticipates" by an amount r the output 
of the master. 

This regime and its stability has been theoretically studied 
in several systems, from the simplest ones described by linear 
differential equations and maps where the mathematical de- 
tails can be fully worked outfl 0Lto the more complicated 
ones such as semiconductor lasers ]6[ operating in the chaotic 
regime. Experimental evidence of anticipating synchroniza- 
tion has been shown in Chua circuits Q and in semiconductor 
lasers with optical feedback HI. 

This same phenomenon has recently been shown to occur 
also when the dynamics, instead of chaotic, is excitable. In 
refs. U the effects of unidirectional delayed coupling be- 
tween two identical excitable systems was studied for both 
the FitzHugh-Nagumo fioll and Hodgkin-Huxley full mod- 
els. It was shown that, when both systems are excited by the 
same noise, and for a certain range of coupling parameters, 
the randomly distributed pulses of the master are preceded by 
those of the slave. This allows for predicting the occurrence 
of excitable pulses in the master. Since many biological sys- 
tems (as neurons and heart cells) exhibit excitable behavior 



and they often operate in feedback regime in a noisy environ- 
ment, the study of the delayed coupling effects in a presence 
of noise is certainly of wide concern. 

The anticipating synchronization regime has been often de- 
scribed as a rather counterintuitive phenomenon because of 
the possibility of the slave system anticipating the unpre- 
dictable evolution of the master |3l 0|. The aim of this 
paper is to provide a simple clear physical mechanism for this 
regime in delayed coupled excitable systems, showing that the 
anticipation of the slave is due to a reduction of its excitability 
threshold induced by the delayed coupling term. As a con- 
sequence, the master and the slave respond to a common ex- 
ternal forcing with different response times. The proposed 
dynamical picture allows us to explain all the general features 
of the phenomenon as well as to determine in a natural way 
the maximum permitted anticipation time. The results are sus- 
tained by numerical integration of the dynamical equations as 
well as by simple analytical calculations. 

A dynamical system commonly used to study excitable be- 
havior is Adler's equation, Jl2ll 



fi — cos x , 



(3) 



This limit cycle develops through an Andronov bifurcation at 
fi c = ±1 II 3H 1411 . where the two fixed points collide and an- 



where x is an angular variable (modulo 2ir) and \i the con- 
trol parameter. For < 1, there are two fixed points at 
x± = ± arccos \i, one being a stable focus and the other 
(x+) an unstable saddle point. If \fi\ > 1, there are no fixed 
points, and the flow consists in an oscillation of the variable x. 
t cycle de\ 

GIB, 

nihilate. For < 1, the system displays excitable behavior: 
if we kick the system out of its stable state with a large enough 
perturbation, the trajectory will return to the initial state (mod- 
ulo 2tt) through an orbit that closely follows the heteroclinic 
connection of the saddle and the node. During this orbit, the 
system is barely sensitive to external perturbations of moder- 
ate amplitude. 

In order to study anticipating synchronization, we consider 
two identical Adler's systems with delayed unidirectional cou- 
pling under the effect of an external perturbation I(t) acting 
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perturbation amplitude po- Note that below the excitability 
threshold, po < 2arccos(/i) (equivalently \x < cos(pn/2)), t r 
does not exist. Forp > cos(po/2) the response time t r is a de- 
creasing function of fj, which approaches zero as \i — > 1. This 
result shows that the response time of an Adler system to an 
above-threshold external perturbation progressively decreases 
as the Andronov bifurcation point = 1) is approached, 
in agreement with the numerical result shown in Fig. 13 (right 
panel). 



FIG. 1 : Time series of the master system x (solid line) and slave sys- 
tem y (dashed line) subjected to white Gaussian noise of zero mean 
and correlations — DS(t — t'), obtained by numerical 

simulation of eqs. I4I5I . Other parameters are: fi = 0.95, K = 0.01 
t — 1. The noise intensity is D — 0.017. 



simultaneously on both systems, 
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When lit) = is zero-mean Gaussian noise, anticipat- 
ing synchronization occurs as shown in Fig. where we plot 
the master and slave outputs for a particular value of K and 
r. Note that the slave system anticipates the firing of a pulse 
in the master by a time interval approximately equal to r. If 
we increase the coupling constant K or the delay time r be- 
yond some values, anticipating synchronization is degraded, 
i.e., the slave system can emit pulses which do not have a cor- 
responding pulse in the master's output, although the reverse 
case never occurs. Upon further increasing K or r, the antic- 
ipation phenomenon disappears. The results are analogous to 
those obtained in 19'] for the FitzHugh-Nagumo model. 

In order to understand the mechanism of the observed phe- 
nomenon, we analyze the behavior of the master alone under 
the effect of a single perturbation I(t) = poS(t — to) acting 
at a certain time to- The effect of this perturbation appears 
only as a discontinuity of the, say, x(t) variable at time to as 
x (to ) = x (to) + Po- The condition for the perturbation to 
be larger than the excitability threshold, is that x(£q") > x+. 
^From now on, we set the initial condition to be in the rest 
state, x(tg) = X-, such the minimum value for the amplitude 
in order to excite a pulse is po > 2 arccos fi and the system 
develops a pulse after a certain response time t r . This time 
can be precisely defined as the time it takes x(t) to reach a 
given reference value, e.g. x r = 7r/2. From Eq. (2) we have 
r*72 — dx — w j 1 j c j 1 yields 
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where b = j^- In Fig. |2](left panel) we plot the response 
time as a function of the parameter fj, for a given value of the 





FIG. 2: Left panel: Response time t r versus /i for the Adler system 
x perturbed by po5(t) with po = 2 from Eq. [6] Right panel: time 
series for x(t) for /i = 0.95 (solid line) and /i = 0.97 (dashed line). 
Both systems have been perturbed at to = 10 by a pulse of constant 
amplitude p = 1.7 and duration At — 0.4. Note that, in agreement 
with the left panel, the system with the larger value of /i pulses before 
the one with the smaller value. 

The fact that the response time decreases with lower ex- 
citability threshold, and that in the coupled system the slave 
can emit pulses that are not followed by a pulse in the mas- 
ter, suggest that the mechanism for anticipation in the master- 
slave configuration is that the slave has a lower excitability 
threshold than the master. This is supported by the following 
qualitative argument: Imagine that at t = to both systems, 
master and slave, are in the rest state x(t^) = y(to) = 
The effect of the perturbation changes both values to x(t$) = 
x(to ) + Po, y{to) = y{to ) + Po- Due to the coupling, 
the slave can be considered to have at this time an effective 
Meff(£o) = M + K[x(t^) - y(t+ - t)] = fi + Kp . Since 
He&(t) > (i also for all times t such that to < t < to + r, the 
excitability threshold of the slave has been reduced and the 
response time decreases. 

To give a more rigorous evidence for this explanation, we 
consider now two coupled systems, Eqs. J4I5L in the presence 
of a single perturbation which we choose to be a pulse of con- 
stant amplitude p and duration At acting at time to in which 
both systems are in the rest state x{t^) = y(t^) = X-. The 
results are reported in figlU 

For a sufficiently large perturbation, the master and the 
slave respond with an excitable spike and the slave pulse an- 
ticipates the master pulse (fig. 3a). For small perturbation 
amplitude no pulses are generated and both systems respond 
proportionally to the applied stimulus (fig 4c). However, an 
intermediate amplitude of the perturbation triggers the emis- 
sion of an excitable pulse by the slave system while the mas- 
ter responds linearly (fig 3b). This confirms a lowering of the 
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I 1 FIG. 4: The ratio between the slave and the master excitability 

r 1 threshold as a function of K for n = 0.05, T2 = 0.2, T3 = 0.35 

; c) J and T4 = 0.5. Considered system have parameter /j, — 0.95. Per- 

: j turbation is applied at time to = 10 with magnitude p = 1.635 and 

10 20 30 40 50 duration At — 0.4. The dashed line corresponds to the constant 

excitability threshold of the master. 



FIG. 3: Response of the master (solid line) and slave (dashed line) 
for three different amplitudes of the singular perturbation of duration 
At = 0.4 at time t = 10: (a) p = 1.7, (b) p = 1.65 and (c) 
p — 1.61. Other parameters are /j, — 0.95, r = 5 and K = 0.01. 

excitability threshold of the slave as compared to the master, 
which is systematically found for all coupling parameters that 
yield anticipating synchronization. In Fig.|4]we plot the ratio 
R between the minimum amplitude of the perturbation which 
generates an excitable pulse in the slave and the minimum am- 
plitude that generates a pulse in the master. As shown in the 
figures, the effect of this particular coupling scheme on the 
slave system is to lower its excitability threshold in such way 
that the difference between the response time of the master 
and the slave to an external perturbation equals approximately 
the delay in the coupling term, r. It is worth noting that when 
K or r tend to zero, not surprisingly the thresholds for the 
slave and the master tend to be equal, while for large values 
of t the difference between the two thresholds is very large. 

Clearly, the same reasoning can be followed if the pertur- 
bation applied to both systems is a white noise process. This 
allows us to explain why the erroneous synchronization events 
correspond to the slave system firing a pulse that is not fol- 
lowed by a pulse in the master: for a particular noise level the 
master response is proportional to the perturbation while the 
slave emits an excitable pulse. By increasing the noise level 
both master and slave emit excitable pulses, each pulse of the 
slave being anticipated respect to that of the master. 

Since, as we have shown, master and slave systems respond 
to external perturbations with different response times, a ques- 
tion which arises is whether it is possible to chose the param- 
eters such that the anticipation time is arbitrarily large, in par- 
ticular, larger than the master response time, r > t r , a result 
that would violate the causality principle. In order to answer 
this question, we plot in Fig. (0 the results of integrating Eqs. 
I4l5l under the effects of a single perturbation for three different 
values of the parameter r. When t < t r or r w t r 105), 



the anticipation time is approximately equal to r. However, 
when t 3> t r , the anticipation time greatly differs from the 
delay time, such that the slave anticipates the master by a time 
interval always lower than t r (jSJ:). This is a reasonable limit 
to the anticipation time: the pulse cannot anticipate the pertur- 
bation which created it. In other words, master and slave are 
both "slaves" of the external perturbation, although the pres- 
ence of the master signal into the coupling term contributes 
to lower the excitability threshold of the slave leading to the 
anticipation phenomenon. 
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FIG. 5: Two coupled systems (master and slave) with a coupling 
parameter K = 0.01 and delay time (a) r = 1, (b) r = 5 and (c) 
r = 50. Both systems have /1 = 0.95 and are perturbed at time 
to — 60 with a pulse of magnitude p = 1.7 and duration At = 0.4. 



In order to assess the generality of our hypothesis, we have 
also considered two delayed coupled FitzHugh-Nagumo sys- 
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tems: 

(2:1,2:2) 

(yi:V2) 



{x 2 + x 1 - -±,e(a-xi)) 



(7) 



fa + yi - f + A >i - fD> £ ( fl - yi)) < 8 > 



In the excitable regime, which occurs when |a| > 1, the 
system possesses a single steady state. As the critical value 
|qj = 1 is approached, the excitability threshold is lowered 
11511 . In this sense, the control parameter a plays the same 
role as the parameter /1 in Adler's equation. In fact, we have 
checked that also in this case the response time of the system 
to an external perturbation decreases as the critical value a c is 
approached (see Fig.|6j. 
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FIG. 6: Time series for the variable x\ of the FitzHugh-Nagumo 
system for a = 1.01 (dashed line) and a = 1.08 (solid line). In 
both cases it is e = 0.09. As indicated by the vertical dotted line, 
the system is perturbed at a time to = 200 by a pulse of amplitude 
p = 0.4 and duration At — 1. Note that the response time decreases 
with increasing a. 

We now consider the unidirectionally delayed coupled sys- 
tem. We find, as in the Adler's system, that the excitability 
threshold for the slave is lower than that of the master, as 
shown in Fig. and that the maximum anticipation time is 
limited by the response time of the master. Finally, we note 
that we have also found exactly the same phenomenology for 
two delayed coupled Hodgkin-Huxley systems. 

The ubiquity of this effect is, in our opinion, an indication 
that the lowering of the excitability threshold of the slave in 
a delayed coupling scheme is a general mechanism for antici- 
pating synchronization in excitable systems. This mechanism 
allows to explain all the observations in the regime of antic- 
ipating synchronization, in particular the erroneous firing of 
pulses in the slave system. In addition, it evidences the causal- 
ity of this phenomenon: the master and slave systems follow 
the applied external perturbations, although the response time 
of the slave system is shorter due to the effects of the cou- 
pling. Moreover, we have shown that the anticipation time is 
limited by the response time of the master system. The rel- 
evance of this type of mechanism for the synchronization of 
coupled chaotic systems is an open question that will be stud- 
ied in the near future. 
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FIG. 7: Response of the master (xi, solid line) and slave (yi, dashed 
line) for two coupled FitzHugh-Nagumo systems with a — 1.01, e = 
0.09, r = 4, K = 0.1, after perturbation at to = 200 by a pulse of 
amplitude p and duration At — 1. For large amplitude, p = 0.4, case 
(a), both systems pulse whereas for the smaller amplitude, p = 0.3, 
there is only pulse in the slave variable. 
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